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Abstract 

We analyze the blowup (finite-time singularity) in inviscid shell models of convective 
turbulence. We show that the blowup exists and its internal structure undergoes a series 
of bifurcations under a change of shell model parameter. Various blowup structures are 
observed and explained, which vary from self-similar to periodic, quasi-periodic and chaotic 
regimes. Though the blowup takes sophisticated forms, its asymptotic small-scale structure 
is independent of initial conditions, i.e., universal. Finally, we discuss implications of the 
obtained results for the open problems of blowup in inviscid flows and for the theory of 
turbulence. 

1 Introduction 

Blowup in a flow of incompressible ideal fluid, i.e., formation of a finite-time singularity in a 
regular solution of flow equations is an important phenomenon, which may be responsible for the 
energy cascade in developed turbulence [191 El HH] • The existence of blowup for fully inviscid 
2D natural convection (Boussinesq equations) and for 3D incompressible Euler equations are the 
outstanding open problems of mathematical fluid dynamics [6l I20|. [2Tj . Dynamical (shell) models 
provide essential simplification of the flow equations, where the Fourier space is substituted by 
a sequence of discrete wave numbers increasing in geometric progression. Despite of simplicity, 
such models preserve many nontrivial features of the original system like, e.g., energy and 
entropy cascades and anomalous scaling of turbulent spectra [4] . Shell models allow very precise 
numerical simulation. This motivates the study of these models with the hope to get a clue for 
explanation of turbulent phenomena. However, turbulent dynamics in shell models is still not 
well understood and explained. 

Mathematical formalism for the Gledzer-Okhitani-Yamada (GOY) [22l [33] and Sabra [26] 
inviscid shell models was developed in [TJ [13], where global existence of weak solutions was 
proved and the criterion for blowup of strong solutions was given. For the dyadic shell model the 
finite time blowup can be proved rigorously [231 19], A constructive method for analysis of blowup 
structure in inviscid shell models was proposed in [15]. This analysis uses renormalized variables, 
where the blowup time is shifted to infinity. Traveling wave solutions of the renormalized 
system lead to the blowup with self-similar structure. Such self-similar blowup's were observed 
numerically in different shell and cascade models of inviscid flows [35[ [3T| [39[ I27j . A similar 
renormalization technique can be introduced for continuous systems [29], see also [16] discussing 
self-similarity of singularities in partial differential equations. 
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In this paper, we consider the shell model of natural convection proposed in [5]. Properties 
of this model in the turbulent regime were studied, e.g., in [37J HI El HD]- in particular, it 
reproduces the Bolgiano-Obukhov scaling law |36[ [25] of turbulent spectra in the inertial range 
(with anomalous corrections). Our goal is to study the blowup in a family of fully inviscid 
models with a parameter controlling the energy transfer between shells. 

First we derive the blowup criterion for the inviscid shell model, which is similar to the Beale- 
Kato-Majda theorem for incompressible Euler equations [3]. Then the renormalization scheme 
is given generalizing the one used in |15[ [29] . Our new contribution is the introduction of the 
Poincare map for the renormalized system. This approach allows rigorous analysis relating 
blowup structures to attractors of the Poincare map. In particular, the self-similar blowup 
structures detected earlier correspond to fixed-point attractors. With a change of parameter of 
the shell model a series of bifurcations is observed, which includes a period doubling cascade 
to chaos and provides a large variety of attractors. We show that all types of attractors (fixed- 
point, periodic, quasi-periodic and chaotic) describe universal asymptotic blowup structures. 
For each type of attractor, we give the detailed description of blowup and determine its scaling 
properties. 

The obtained results give some insight to the open problems of blowup in continuous inviscid 
flow models. We see that the blowup can have the universal asymptotic structure, which does 
not depend on initial conditions up to system symmetries. However, this structure does not have 
to be self-similar and a number of different scenarios including periodic and chaotic regimes are 
possible. Thus, it may be very hard to detect such universality, unless special renormalization 
tools are applied like the one proposed in this work. 

The paper is organized as follows. In Section 2 we introduce the model and give the blowup 
criterion. Section 3 describes the renormalization scheme of [15], which is demonstrated on the 
isothermal shell model in Section 4. The Poincare map for renormalized dynamics is defined in 
Section 5 and the corresponding bifurcation diagram is studied in Section 6. Sections 6-9 analyze 
blowup structures corresponding to fixed-point, periodic, quasi-periodic and chaotic attractors 
of the Poincare map. In Section 10 we summarize the results and discuss their relation to the 
open problems of blowup in continuous models and anomalous spectra of convective turbulence. 



2 Blowup in inviscid models of natural convection 

In the Boussinesq approximation |24] , buoyancy-driven flows of unit density are governed by the 
equations 

0u d0 

— = -u- Vu- X7p + vAu + ag9e z , — = -u ■ V9 + kA9, V • u = 0. (1) 
at at 

Here u is the velocity, p is the pressure, and 9 = T— Tq is the deviation of temperature T from the 
mean value To. The constant parameters k, v, a, and g are the thermal diffusivity, the kinetic 
viscosity, the thermal expansion coefficient, and the gravitational acceleration, respectively; the 
unit vector e z specifies the vertical direction. In the inviscid limit, the system conserves the 
entropy S = f ^9 2 dV and the total energy E = J(i|u| 2 + agz9)dV, where z is the vertical 
coordinate. 

Our work is motivated by the problem of blowup of classical solutions for the fully inviscid 
system, i.e., formation of a finite-time singularity in smooth solutions when v = k = 0. In the 
isothermal case, 9 = 0, inviscid equations ([T]) reduce to the incompressible Euler equations. It 
is well-known that the Euler equations for ideal incompressible fluid in 2D possess a global in 
time, unique, regular solution, see, e.g., [30]. Global regularity of the 3D incompressible Euler 
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equations is an open problem [6j[20l[2T]. The Beale-Kato-Majda theorem [3] states that, if the 
initially smooth solution of the Euler equations cannot be continued beyond the time t c and t c 
is the first such time, then 

rt c 

||w(-,t)||oodt = °o, (2) 



where ||u>(-, f)||oo denotes the L°° norm of the vorticity u = V x v at time t. This condition 
implies that the maximum vorticity must grow at least as (t c — near the blowup. In the 
general case 9^0, the existence of blowup in smooth solutions of the fully inviscid Boussinesq 
system ([I]) is an open problem both for 2D and 3D flows [6]. The blowup criterion for 2D flows 
includes ([2]) and an extra condition for the temperature gradient |40t [7]. 

Dynamical models, also called shell models, represent an essential simplification of system ([I]) 
constructed in the discretized Fourier space k n = koh n with h > 1. Here k n with n = 0, 1, 2, . . . 
is the wavenumber, and the corresponding effective velocity and temperature are given by real 
quantities u n and 9 n . We will consider the shell model proposed in [5], which has the form 



~dT 

d6n 
dt 



Ak n {u n _i — hu n u n+ i) + Bk n {u n u n _i - hu 2 n+ i) - vk 2 n u n + ag6 n + f n , (3) 
A& n (it n _i0 n _i - hu n 9 n+ x) + Bk n (u n e n -x - hu n+1 9 n+1 ) - nk^9 n + g n , (4) 



where A, B, A, and B are real parameters and n_i = = 0. Equations ([3]) and have 
the structure similar to the Fourier transformed Boussinesq equations (pQ). Nonlinear quadratic 
terms describe the shells interaction and involve the nth and its neighboring shells. The addi- 
tional forcing terms f n and g n model the influence of boundaries (e.g., due to heating from below) 
and are usually restricted to the first shells. Despite of simplicity, the shell models demonstrate 
several highly nontrivial properties of the original continuous system ([T]), for example, the tur- 
bulent regime with entropy cascade and the anomalous statistics close to the Bolgiano-Obukhov 
scaling law in the inertial range [5j [37] . 

We assume the traditional choice of parameters, h = 2, ko = ag = B = A = B = 1, and 
consider A = e. The parameter e controls the nonlinear process of energy transfer between 
shells, and the value e = 0.01 was chosen in earlier works, e.g., [10] . We will study the fully 
inviscid equations ([3]), (J1J) with no forcing, which reduce to 

^ = k n [s(Un-i ~ hu n u n+1 ) + n n ii n „i - hu 2 n+l ] + 9 n , (5) 



dt 

dBn 

dt 



kn (u n —\9 n —\ - hu Tl 9 n+ i + u n 9 n _i - hu n+ i9 n+ i) . (6) 



This system conserves the entropy S = 2®n (here and below the sums are taken over all shells 
n). If 9 n = 0, the energy E = ^ ^u n is conserved. However, there is no analog for the total 
energy in the case 9 n ^ for this shell model. 
Following [13], we introduce the norms 

ll^" 1 := (y2 k n u l) i IMIi.oo := sup&„|u n |, (7) 
and consider smooth solutions with finite norms 

|H|i + ||0||i < oo. (8) 
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Local existence of such solutions can be proved using the Picard theorem in the same way as for 
the inviscid Sabra shell model [13] with obvious modifications due to different nonlinear terms. 
The blowup at t = t c implies that 

sup (||n||i + ||#||i) = oo. (9) 

o<t<t c 

Note that the blowup does not imply a singularity for a particular shell variable u n or 9 n , 
since they represent Fourier components of the flow; only the integral contribution (the norm) 
explodes. The following theorem provides the blowup criterion, which is similar to (|2|) and to 
the one derived in [T3] for the Sabra shell model. 

Theorem 1 Let u n (t) and 8 n (t) , < t < t c be a smooth solution of the inviscid shell model fSP, 
satisfying condition where t c is its maximal time of existence. Then, either t c = oo or 



/ |M|l,ooCtt 
J 



co. (10) 



Proof: Condition (|10p implies that ||it[|i oo is unbounded for < t < t c and, hence, condition 
([9]) is satisfied. Therefore, (|10p is a sufficient condition for the blowup. 
Using (jSJ), ([6]) and the norm defined in (J7J), we derive 

\dt + "^"^ = ^ Unk ^ [ £ ( U n-± ~ hu nU n +l) + «n«n-l ~ ^n+l] + ^ k^U n 6 n 

+ ^2 On^n ( u n~lS n ~l ~ hu n 9 n+ i + U n Q n -\ — hu n+ i9 n+ i) . 

(11) 

Using the condition k n \u n \ < \\u\\i t00 and applying Cauchy-Schwarz inequality in (fTTI) yields 

j i\H\ + < C||«||i )00 ||u||? + 2||«||i||e||i + C||u||i )0 o||e||? < (1 + C\\u\\ lj00 ) (\\u\\\ + \\e\\\) 

(12) 

valid for a sufficiently large positive constant C. Using the Gronwall's inequality we get 
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i + ||^||i) t=tc < (||n|| 2 + ||^|| 2 ) 4=0 exp^ c + C^ C ||n|| liOO d^ . (13) 

This proves that condition (|10p is necessary for the blowup. □ 

Condition (|10p implies that the maximum shell "vorticity" must grow at least as k n u n ~ 
(t c — t)^ 1 as t —¥ t~ . It is interesting that no extra condition for the shell temperatures is 
necessary. We remark that the statements of this section (on the local existence and blowup) 
remain valid for the sequences with n £ Z, i.e., without the left end at n = 0. Also, the norm 

1 /2 

||n||i can be replaced by \\u\\d = (Y^ fc 2d u 2 ) for the sequences with n > 0, but not for n £ Z. 



3 Dombre-Gilson renormalization scheme 

In this section we introduce the renormalization scheme analogous to the one suggested by 
Dombre and Gilson |15j for the mixed Obukhov-Novikov (ON) model |32|. 114], The aim of this 
renormalization is to obtain a system, where the blowup time t c is moved to infinity, so that 
standard dynamical system methods can be applied. 
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The new time variable r is denned implicitly as 

t = jT exp (- jT A{r")dT''\ dr', (14) 

where the function A(t) will be specified later. The scaled speed v n and temperature $ n are 
given by 

V n = exp (- j A(T')dT^j k n U n , ti n = exp ^-2 J A{T')dT^j k n 9 n . (15) 

Differentiating expressions (|15|) with respect to r and using (JSJ), (0) and (|14|) yields the system 
^ = P n - Av n , P n = e(/i 2 ^_ 1 - v n v n+ i) + hv n v n -i - h~ x v 2 n + x + i9 n , (16) 

—j^- = Q n - 2Av x n , Q n = /i 2 W n _ilV-l - t>„#n+l + /l^n-l ~ h~ . (17) 

It is easy to see that 
d 



dr 



£(«n + O = + l^nQn) ~ 4A + (18) 



Let us choose 



Then the sum 



Mr) = (v n P n + h n Qn^j /E K + <£) • 



(19) 



£>*+<) = c (20) 

is conserved. At r = expressions (fT5]) and 0, (jSJ) yield ^ v n < (^v^) 2 = \\u\\f < oo and 
^2$ n = \\0 ||f < co. Hence, the sum in (|20p is finite, c < oo. In the next lemma we show that 
the function in (|19p is well defined for any nontrivial solution (we will not consider the trivial 
solution v n = 0, i? n = from now on). 



Lemma 1 Nontrivial solution v n {r), $ n ( T ) of system M6\) . |-/7| ) exists globally for < r < oo 
and it is related by {L$ , / T?5j) to i/ie solution u n (t), 6 n (t) of system |5|), |j|) with t < t c , where 
t c is the blowup time from Theorem [I] 

Proof: Because of relations ()14|) . f)15|) between the solutions, we need to show only that A(r) in 
(fT9|) is well defined and that any finite r > corresponds to £ < t c . Due to (|201) . we have 

|«n| < c 1/4 . (21) 



The denominator in (119]) is constant. We substitute P n , Q n from (fTBjl . (1171) into (119]) and consider 
the terms in the nominator. The sum of first terms can be bounded as 

^ 2 £«-i] < k!^c 1/4 E«_i < \e\h 2 c^J2vi < (22) 

where we used (121 p . the Cauchy-Schwarz inequality and (I20p . Similar bounds can be found for 
the other terms. It follows that |A(r)| is bounded for all r > 0. 

Using (I2ip and finiteness of the integral in (115p . one can show that the sequence \k n u n (t)\ is 
bounded for any fixed t corresponding to < r < oo. Hence, |[u||i >00 < oo and, by Theorem (TJ 
we have t < t c . □ 
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Equations (fT6|h (jTTJ) with A from (fl~9j) possess three types of symmetries 

(i) r !->■ r/a, t> n H> au n , i? n (->• a 2 tf n ; (23) 
(ii) r i-> r - r ; (24) 
(iii) f n l-> Vn+l, i? n H> (25) 

The last symmetry does not hold for the equations of the boundary shell n = 0, but this 
symmetry becomes exact when we pass to the infinite lattice n € Z in Section [SJ 

Lemma 2 XTie symmetries A23\) -[25 \) are equivalent to the following symmetries for the original 
shell model (0): 

(i) t*-¥t/a, u n \-¥ au n , 9 n ^a 2 9 n ; (26) 
(ii) th^t/a — to, u n h->- au n , 6 n ^a 2 6 n , (27) 
(Hi) u n ^hu n+1 , 9 n h^h9 n+1 . (28) 

In the case (ii) the values of to and a are uniquely determined by tq. 

Proof: Symmetries (|26|) - (|28|) describe time scaling, time shift and space scaling (recall that 
k n+ i = hk n ), and can be related to ([23|) - (f2"5|) using (fT4"J) and (fT5j) . We will consider in detail the 
most difficult case (ii). The symmetries (i) and (Hi) are studied similarly. 
According to (|24p . the new variables (denoted by a hat) are 

t = t To, v n (?) = v n (r), d n (?) = i? n (r). (29) 

It follows from 032) that 



A(t) = A(t)=A(t + t ). (30) 
For the time t corresponding to t, expression (114p yields 



(31) 



/■T / /V ^ \ i-t—to ( pt' \ 

t = y exp f - J A(T")dr" JdT> = J ex P ( " / A ( T " + T °) dr " j dr ' 

= fexpf- f A(OdAdt\ 

where we used r = r — tq from (|29p , A from (|30p and made the changes of variables £' = t' + to 
and £" = r" + tq. Expression (|3ip can be written as 



t = t/a- 1 , (32) 

where 

a = exp (- ° A(T")dr'^j , to = £° exp ^- jf A(t")cZt"J cfr' (33) 
and i is given by {Ml). For the shell speeds, the first expression in (fT5"j) with (f2Uj) and (f3"Uj) yield 

u n (i) = exp ( jf l(r')dr' J A^^nCf) = exp ° A(t' + t )(Zt'J fc~ V(r). (34) 
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Figure 1: Blowup in isothermal (mixed ON) shell model: (a) traveling wave in renormalized 
variables v n (r), (b) universal function U(t), (c) numerical simulation of shell velocities u n (t) 
(black curves) and their asymptotic self-similar form (bold light-red curves). The blowup time 
t c = 2.07 corresponds to r — > oo. 

The change of variable £' = r' + To in the integral yields 

u n (t) = exp ^ A(T')dT^j k^Vnir) = au n (t) (35) 

with a from (|33p . Similarly, one derives n (t) = a 2 9 n (t) using the second expression in (115p . 
which accomplishes the proof of (|27p . □ 

4 Universal structure of blowup 

As a demonstration of the renormalization method, let us consider the isothermal (mixed ON) 
shell model with 6 n = $ n = and e = 0.5. Numerical solution of equations (|16j) is shown 
in Fig. []Ji. Here we used 100 shells and the initial condition with only the first shell speed 
perturbed, i>o(0) > 0. For large times, the asymptotic solution has the form 

Vu{t) = aV(n — ar), (36) 

which was first observed in [15]. Expression (|36p defines a solitary wave traveling to larger shell 
numbers n with constant speed a. The function V(£) vanishes in the limits £ — )• ±oo, and the 
constant a > is arbitrary reflecting the scaling symmetry (|23p . 

The theory explaining asymptotic solution ([36]) will be given in the next sections using 
the dynamical system approach. This solution yields the self-similar expression for original 
variables u n (t). In the next theorem we provide such self-similar expression and its derivation 
(in modified form) given by Dombre and Gilson; we will need these results later for the study of 
nonzero shell temperatures. Note that analogous self-similar blowup was observed numerically 
in [35\ [3"H [39| [27] for different dynamical models of hydrodynamic flows. 

Theorem 2 ( |15j ) Taking a = 1 in $36\) we define the real quantity 



1 



" = Whk A{T)iT (37) 
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and the function 

U(t-t c ) = exp(J A{T')dr^jV{-T), (38) 

where r is related to t implicitly by {1$ and A is given by U9\) . If y > 0, then the solution u n {t) 
corresponding to A36\) has the form (for arbitrary a > 0) 

u n (t) = aky- 1 U(aky(t-t c )), (39) 

where 

t c = jf exp ^- ^(r'Odr'^ dr' < oo (40) 

zs i/ie blowup time. 

Proof: Let us consider the time-scaling symmetry ([23j). By Lemma [H the corresponding quan- 
tities t, t c and u n scale as 

i-i c i-4 (t-t c )/a, u n >-^au n . (41) 



This shows that dependence on a in (|36p and (|39p is the result of the same symmetry. Hence, 
we can take a = 1 in the rest of the proof. 

The function A{t) given by (|19p and (|36p is periodic with period 1/a = 1. Then, it follows 
from (I37p that the inequality 

^(r)dr > C + ry log /i (42) 



holds for some constant C. Using this estimate, it is straightforward to show that the integral 
in (|40|) converges for y > 0. 

Due to periodicity of A(t) and ([37|) we have 

exp / A{r')dT' = fcv (43) 



for any r, where fc n = /i n . Let us consider the time t' corresponding to r + n. Using (|14[) and 
(SOl), we obtain 



/oo / /-t' \ poo ( pr+n \ 

exp f - y A(r")^r" j dT ' = J exp ( " / A(r")dT" j (if 

= exp ( - jT A(r")dr" J = k?(t e - t), 



(44) 



where we changed the integration variable t' = t + n and then used (|43p . Similarly, we express 
the shell speed u n from (fT5|) and use (03]), (|36|) and (|38|) as 

-u„(t') = A;" 1 exp ^jf + A(rOdr^ w n (r + n) = k v ~ x exp A(r')rfr'^ U(-r) = k y n ^ l U{t-t c ). 

(45) 

Substituting t — t c = knit' — t c ) expressed from (j4"4"|) into (j4"5|) and dropping the prime, we obtain 
relation (|39|) with a = 1. 

Finally, let us show that expression (|39|) describes the blowup at t c . For this purpose we 
consider the times t n = t c — bkn V with arbitrary b > such that U{—b) ^ 0. Then expression 
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(f39|) with a = 1 yields k n u n (t n ) = knU(-b). For y > we have fe n |« n (t n )| — ► oo as n — > oo and, 
hence, the norm ||u||i —> oo. Since t n — > t c , this means the blowup at t c . □ 

Since (|36p is an asymptotic solution for the renormalized system, expression (|39p represents 
the asymptotic form of blowup for large n and t — >■ t~. In our numerical example, computations 
yield y = 0.8557 and the function is shown in Fig. [lb. Figure [lb compares the asymptotic 
self-similar relation (|39p with the numerical shell speeds u n (t) demonstrating convergence for 
large n. The scaling exponent y and the function U (t) in expression f|39|) do not depend on initial 
conditions, which was confirmed numerically and will be explained theoretically in Section [71 
Therefore, the blowup has the universal self-similar asymptotic form. 

5 Poincare map 

Asymptotic traveling wave solutions (|36|) do not exist for all values of parameter e. Instead, 
solutions of system (|16p . (|17p appear in the form of periodically or chaotically pulsating waves. 
We propose a general method for analysis of such solutions following the dynamical system 
approach and using the Poincare map introduced below. 

Let us consider the infinite-dimensional space W of renormalized flow variables 

W = (. . . ,V n -i,V n ,V n+1 , . . . ,1? n -l, fin^n+l, ■■■) EW (46) 

with the £ 2 norm ||w|| 2 = ^2(v 2 + i? 2 ). According to (|15h and d7j, we have ||w|| 2 = c u ||«|| 2 + 
c e ||$ Hi i where the norms and the coefficients c u , eg are finite for any r > by Lemma [TJ Hence, 
w(r) € W for r > 0. 

The blowup phenomenon is described by dynamics for large shell numbers n. Hence, we 
can disregard the left end (n = 0) of the shell sequence in the blowup analysis, and consider 
equations (|16p . ()17p for the sequences with n G Z. It is easy to check that all the statements 
made above about the shell model solutions remain valid in the case n € Z. 

We define the real number 

n„(r) = £ n + ^) («* + , (47) 

which estimates the shell reached by the flow at time r. Recall that the denominator in (|47p is 
constant, see ()20p . Let us consider w = w(0) as the initial condition for r = 0. Then we define 
the vector w' = w(ti) as the solution at the first time moment n > when the position n„ 
increases by one, i.e., 

n„(ri) = ra„(0) + 1. (48) 
The corresponding transfer operator is denoted by 71, so that 

Tiw = w', v' n = v n (ri), K = Mn)- (49) 

The operator 7i is defined for given w if condition (|48j) is satisfied for some ri > 0. Numerical 
computations presented below suggest that this operator is defined for any nontrivial vector 
w E W, i.e., for any nontrivial solution we have n v — > oo as r — > oo. 
Also, we define the shift operator S : W h-> VF acting as 

5w = w', ^ = w n+ i, i?^ = i? n+ i. (50) 

This operator represents the renormalization in physical space, because n + 1 corresponds to the 
wave numbers k n+ \ = hk n . 
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Figure 2: Numerical simulations for the isothermal shell model with e = 0.5. (a) Discrete 
dynamics induced by the transfer operator T\- Shown are the shell speeds v n at times r m 
corresponding to the wave positions n v (T m ) = m for m = 0, 1, 2 . . . (b) After the shift S m to the 
left, these solutions represent iterations of the Poincare map V m and converge to the fixed-point 
attractor (bold red curve) as m — > oo. 



Now we can define the Poincare map as V = 571, i-e., 

•pw = w', v' n = v n+ i(T 1 ), ■&' n = 'd n+1 (ri). (51) 

Therefore, the Poincare map V represents the time-transfer operator combined with the shift of 
shell numbers. Due to this shift, the Poincare map describes the flow dynamics observed in the 
"moving frame" along the logarithmic axis n = \og h k n in the Fourier space. 

In the isothermal shell model of previous section, solution (|36[) corresponds to a fixed-point 
attractor of the Poincare map, see Fig. [2j We will see that the case of nonzero shell temperatures 
yields variety of attractors depending on the parameter e. Note that the approach proposed 
in this section is similar in spirit to the renormalization group method in partial differential 
equations, see, e.g., [HJ [2]. 

6 Bifurcation diagram 

We studied numerically the limiting behavior of a discrete dynamical system governed by the 
Poincare map (|5ip . For this purpose, we used 100 shells in system (|16p . f)lT[) . This gives a very 
good approximation, because the solution decays rapidly on both sides and occupies effectively 
about 10 shells, see Fig.[2j Initial conditions with a few nonzero components centered at the shell 
n v = 70 were chosen. In order to remove the transient behavior, the first 4000 iterations were 
skipped. Figure [3] shows the bifurcation diagram represented by the speeds v n and temperatures 
d n of the shell n v = 70. 

Fixed-point (period-1) attractors exist for e < —1.255. At e = —1.255, the period-doubling 
bifurcation occurs leading to a period-2 attractor. The next period-doubling bifurcation corre- 
sponds to e = —0.749, giving rise to the period-4 attractor. Further increase of e leads to the 
infinite period-doubling cascade, which follows the Feigenbaum scenario [18] and ends with the 
chaotic behavior for e > —0.7. For larger e a number of windows with periodic solutions are 
found, which separate the regions with chaotic or quasi-periodic attractors. The largest periodic 
windows are indicated in Fig. [3l At e = 0.015 a subcritical bifurcation occurs (see Section [7]) 
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Figure 3: Bifurcation diagram shows (a) the renormalized shell speeds v n and (b) the renormal- 
ized shell temperatures ~d n on the attractor of the Poincare map for n = n v . The fixed-point 
attractors F + , F_ and F c are distinguished, see Section [TJ The largest periodic windows P^, are 
indicated, which correspond to period-ra attractors of the Poincare map. 

leading to a fixed-point attractor for e > 0.015. Numerical results suggest that there is a unique 
attractor for each value of e (up to system symmetries), but multiple attractors may appear in 
general. 

Examples of periodic, quasi-periodic and chaotic attractors on the plane (w n ,$n) are shown 
in Fig. [U In the next sections we describe how each type of attractor determines the specific 
blowup structure. 

7 Self-similar blowup solutions 

First, let us consider a fixed-point (period-1) attractor w of the Poincare map, i.e., 

Tw = w. (52) 

In this case t\ = 1/a, where 1/a is the time period determined by the Poincare map. Using 
definition ([5"T]) for the shell speeds in ([52]) . we have 

v n+ i{r + 1/a) = v n (r). (53) 

It is easy to see that a general sequence v n {r) satisfying condition ([53]) has the form of a traveling 
wave (]36j) . The corresponding solution for the original variables u n {t) is described by Theorem[2l 
Figure [5h. shows the scaling exponent y from (|37p as a function of e computed numerically (the 
values of y corresponding to other types of attractors are discussed later). Since y > 0, the 
solution blows up at finite time t c given by (|40p . 

All shell temperatures vanish, = 0, for the fixed-point attractor, as one can see from 
Fig. [3b. In order to capture the limiting behavior of the original variables 9 n (t) near the 
blowup, we must consider small perturbations ow(r) near the fixed-point w. Following the 
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Figure 4: Periodic, quasi-periodic and chaotic attractors on the plane (f n ,i?n) for n v = 70. 
These attractors correspond to e = —0.5, e = —0.35, and e = 0.01, respectively. 
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> - 
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Figure 5: (a) Scaling exponents y (black) and z (red) as functions of the shell model parameter 
e. The exponent z corresponds to fixed-point attractors. (b) The largest magnitude Floquet 
multiplier p of the fixed-point attractor depending on e. 
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classical perturbation theory (see, e.g., [M]); we consider linearization of the Poincare map near 
the fixed point, 

V(w + Sw) taw + dV Sw, (54) 

where dV is the linear part (Jacobian matrix) of V . Let p be the Floquet multiplier and w be 
the corresponding eigenvector satisfying the eigenvalue problem 

dV w = pw. (55) 

The multipliers are complex numbers in general and satisfy the condition \p\ < 1 because w is 
an attractor. The asymptotic behavior of small perturbations 5w(r) for large times is given by 

5w{T m ) « (dV) m 5w{0) « p m w, (56) 

where r m = m/a and p is the Floquet multiplier with maximum absolute value \p\. A linear 
combination must be taken in f|56H if there are several p with equal maximum absolute values. 

Since # n = for the fixed-point, the perturbed solution is # n = 5$ n and we will drop 5 
below. The relations (15ip are linear and, hence, are valid for the linearized operator dV as well. 
Therefore, similarly to ([52]) and ([53]) . the eigenvalue equation ([55]) for the shell temperatures 
yields 

tf n+ i(r + l/a) = P K(t). (57) 
A general sequence $ n (r) satisfying (j57p can be written in the form similar to (|36p as 

Mr)=bp n Q(n-aT), (58) 

where 6 is an arbitrary constant. 

The original shell variables 9 n (t) are expressed from (I58p using the formulas (j!4p , (I15p , where 
the quantities r and A correspond to the unperturbed (fixed-point) solution. The derivations 
analogous to those made in the proof of Theorem [2] for shell speeds (|36|) yield 

9 n (t) = bp^-^akUt - t c )). (59) 

Assuming a = 1, the function is determined by the expression similar to (|38|) as 



Q(t-t c ) = exp ^2^ ^(r')dr'J Q(-r). (60) 

The Floquet solutions and multipliers can be found numerically by standard methods, see, 
e.g., [31]. The multiplier p with maximum absolute value depending on parameter e is shown in 
Fig. [5b. The multiplier is positive, negative or complex in three intervals denoted, respectively, 
by F+ = [0.5, oo), F_ = (-oo, -1.255] and F c = (0.017, 0.5), see also Fig. [3] In the tiny 
interval 0.015 < e < 0.017 (not shown in the figures), a real negative multiplier gets larger in 
magnitude than the complex ones. At e = 0.015 a bifurcation occurs {p = —1). Since the 
fixed-point solution becomes unstable and no small-amplitude attractor appears for e < 0.015, 
this bifurcation is subcritical, see, e.g., [38] . 

We define 

\p\ n = k- z , z = - log h \p\ > 0, V3 = argp. (61) 

The scaling exponent z depending on e is shown in Fig. [5^. For real positive p, expression (]59|) 
takes the form 

p>0: 6 n {t) = bk 2 y- z - 1 Q(akl{t-t c )). (62) 
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Figure 6: Thin black curves show numerical solutions 9 n (t) for (a) e = 0.7, (b) e = —1.5, (c) 
e = 0.3. These solutions are governed by the Floquet states with positive, negative and complex 
multipliers, respectively. Bold light-red curves represent the universal asymptotic form of the 
blowup given by (|62p -(|64 p . In order to demonstrate convergence for large n, several initial shells 
are skipped. 



In the case of real negative p, we have 

p < : 9 n (t) = b(-l) n k%>- z - l <d{aky n (t - t c )). (63) 

When p is complex, the function = ©i + i©2 an d factor b = b\ + ib 2 are also complex, and 
the real solution is obtained by taking real part of (159p . This yields 

complex p : 9 n {t) = k^'^Re + ib 2 )e in * '(©i(£) + i© 2 (£))] > £ = ak*(t - t c ). (64) 

Figure [6] presents numerical results for e = 0.7, —1.5 and 0.3 corresponding to positive, negative 
and complex Floquet multipliers p. The numerical solutions 9 n (t) are shown (black curves) 
together with their asymptotic form (bold light-red curves) determined by expressions (|62|) — 

§3). 

Since the fixed-point ([52]) is an attractor, asymptotic solutions (|36l) and (i58]l are determined 
uniquely up to system symmetries ([23]) - ([25|) . For the original variables u n (t) and 6 n (t), these 
symmetries have the form (I26p -( i28l) . One can check that these symmetries reduce to a change 
of the coefficients a and b in ([39]) and ([62]) - ([M|) only, leaving the scaling exponents y, z and the 
functions U, unaltered. We confirmed this fact numerically as well. 

We conclude that the fixed-point attractor of the Poincare map V determines the asymptotic 
self-similar form of blowup. This form is described by expression (|39p for shell speeds and by 
(|62p ~ (|64p for shell temperatures, where a > and b £ M. are arbitrary amplitudes. The scaling 
exponents y and z, as well as the functions U and do not depend on initial conditions for a 
given shell model. Therefore, the asymptotic self-similar form of blowup is universal. 

8 Blowup with periodic structure 

Let us consider a periodic attractor of the Poincare map given by the vector w such that 

V p w = w, (65) 
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where the period p > 1 is a positive integer. If p/a is the time period of V p , then the times r n 
corresponding to n iterations of the Poincare map satisfy the conditions 

T n = Tj+pN/a, n = j + P N, j = 0,...,p-l, N = 0,1,2,... (66) 

Using (f66|) and ([5T]). periodicity condition ([65]) can be written as 

v n + P {r + p/a) = v n (r), tf n+p (r + p/a) = i? n (r). (67) 

A general form of the functions satisfying (|67p is 

« n (r) = aVj(n - ar), i? n (r) = a 2 QAn - ar), j = 0, ...,p-l, (68) 



where n = j + piV and the factors a and a 2 are used to reflect the symmetry (|23p . For each j 
solution (|68p represents a wave moving with speed a along the logarithmic axis n = \og h k n . 



y = —^— f A(r)dr (69) 



Theorem 3 Taking a = 1 in I168\) we define the real quantity 

1 

plog h 



o 



and the functions 



kj v exp A^dr'JVjij-T), 



Uj(t - t c ) = k~ y exp / A(r')dT' Vj(j - r), (70) 



&j(t-t c ) = kj 2y exp(2^ A^dAQjij-r), (71) 

where t is related to t by (T4\ ) and A is given by [TS\). If y > 0, then the solution u n (t), 9 n (t) 
corresponding to $68\) has the form (for arbitrary a > 0) 

u n (t) = akl^U^akl^t - t c )), Mt) = « 2 ^" 1 © J «- j (t " *c)), (72) 
and describes the blowup at t c given by ( f^0| ). 

Proof: The coefficient a in (|68p and (|72p reflects the time-scaling symmetries (|23|) and ()26p . 
Hence, we can take a = 1 in the rest of the proof. We will derive expressions (|70p and (|72|) for 
the shell speeds; the shell temperatures in (|7ip and (|72p are considered similarly. 

The function A(r) given by (fT9|) and (f68l) is periodic with period p. Due to (j69j) we have 

/ f-T+pN \ 

exp^jf A(r')dr'j = k y pN (73) 

for any r, where fc p 7v = h pN . Let us consider the time t' corresponding to r + pN. Using (|14p 
and (j40|) . the derivation analogous to (f44j) yields 

t c -t' = k^(t c -t). (74) 

From (I68p we find v n (r + pA^) = Vj(n — r — pN) = Vj{j — r) for n = j + pN. Then we express 
the shell speeds u n from (fT5|) and use ([73]) . (f70|) as 



/ rr+pN \ 

tin(t') = k~ x exp ( / A(r')dr' w„(r+pJV) 

= k-^expU A(r')dr'y 3 (j -r) = kl~ l U 3 {t-t c ). 
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Figure 7: Asymptotic form of renormalized shell speeds v u (t) and temperatures $ n (r) for the 
blowup with periodic structure (j68[) in the shell model with e = —0.5. This solution corresponds 
to the P4 window in Fig. [3l 

Substituting t — t c = k y n _j(i! — t c ) expressed from (jJ3J) into (jT5"j) and dropping the prime, we 
obtain the first formula in (j72|) with a = 1. 

It follows from periodicity of A[r) and (|69p that the inequality (|42p holds for some constant 
C, which implies convergence of the integral in (|40p and the blowup of solution (|72p at i c for 
y > 0, just as in the proof of Theorem [2j □ 

Since $ n 7^ in the periodic attractor, the speed and temperature variables are strongly 
coupled in the blowup, in contrast to the case of the fixed-point attractor in Section [71 As 
a consequence, there is a single scaling exponent y in self-similar expressions of Theorem [3l 
Periodic solution (|65p is an attractor and, hence, expressions (|68p represent the asymptotic form 
of solution in renormalized variables for large t, Fig. [7J The functions Vj and Qj in (|68p are 
determined uniquely up to system symmetries. The symmetry (j23|) is taken into account in (|68p . 
and the symmetries ()24p . ()25p reduce to the shift of time r and the cyclic permutations of the 
indexes j i-> j + 1 (mod p) . 

We conclude that the blowup with periodic structure has asymptotic form (172p for large n and 
t — ?• i~, where the scaling exponent y and the functions Uj, Qj are universal, i.e., independent 
of initial conditions (up to cyclic permutations of indexes j = 0, . . . ,p — 1). The value of y 
depending on the parameter e is shown in Fig. [5^. Figure [8] presents a numerical example of the 
blowup. 

9 Blowup with quasi-periodic and chaotic structure 

According to numerical simulations in Figs. [3] and U] quasi-periodic and chaotic attractors exist 
in the interval —0.7 < e < 0.015, outside the periodic windows. In particular, the value e = 0.01 
used in [5j [TO] corresponds to the chaotic attractor. Figure Eh shows the times r n and the 
integrals f Q n A(r)d,T corresponding to n iterations of the Poincare map for e = 0.01. These 
quantities grow linearly with n up to small chaotic oscillations. 
We introduce the quantities 



- = lim ^ > 




y = 



a log h 



(76) 
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Figure 8: Shell velocities u n (t) and temperatures 9 n (t) for the shell model with g = —0.5. The 
shells n = 24, . . . , 34 are shown, which correspond to the solution displayed in Fig. [7j The 
solution blows up as t — > t~ and has the period-4 structure visible in the figure. 




Figure 9: (a) The times r n and integrals I n = f n A(t)cIt corresponding to n iterations of 
the Poincare map V for g = 0.01 (chaotic attractor). The dashed lines indicate the slopes 
1/a and (A) /a = ylogh. (b) Power law scaling of shell speeds u n (t) and temperatures 9 n (t) 
corresponding to the times when n v = n. The dashed lines indicate the slopes y — 1 and 2y — 1 
given by ([78|) and ([79]). 
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where 1/a represents the mean time step (r n — r n _i) of the Poincare map and (^4) is the mean 
value of the function A{r) on the attractor. Note that an arbitrary value a > can be obtained 
after the time-scaling symmetry transformation (|23p . but this transformation does not alter the 
value of y. Also, for fixed-point and periodic attractors, the last expression in (|76p reduces to 
(1371) and p9l . 

Since A(t) = (A) + 5A(t), where 5 A oscillates near zero mean value, the inequality (|42p 
written for a = 1 holds, see Fig. [9k.. As a result, the integral in (flO|) converges for y > 
providing the finite blowup time t c . 

Figure [TOa shows the speeds v n (T), whose dynamics can be viewed as a chaotic wave moving 
toward large shell numbers n. Since one iteration of the Poincare map corresponds to the increase 
of wave position by one, the mean wave speed equals a. The wave amplitude in renormalized 
variables v n and # n does not change, see (|20p . Using (|76p we estimate 

exp (J ^ A(t')cIt'^ = e {A)Tn exp QT" SA^dr^ ~ e^ T " ~ e (A)n/a ~ fc» (77) 

Therefore, when the wave reaches the shell n, the shell speeds and temperatures can be estimated 
by order of magnitude using (fl"5j) as 

u n = k^expU^ A(T')dAv n ~ky-\ (78) 

9 n = k^expM^ AiryAtin-k'*- 1 . (79) 

The corresponding time interval to the blowup is found using (|14p . ()40p . (I76p and (]77p as 




where we changed the integration variables as r' = £' + r n and r" = £" + r n ; note that the last 
integral factor represents an oscillating quantity with finite mean value. 

Scaling laws ()78|) — ()80|) are the same as for the blowup with periodic structure, see (|72p . They 
are confirmed numerically in Fig. [9)3, which presents the logarithmic plots of the amplitudes 
|u n (t)| and |0 n (i)| & t the times t when the wave position n v is equal to n. The scaling exponent 
y computed numerically is shown in Fig. [5^. Since y > 0, the estimate (I78p shows that k n u n 
becomes infinitely large as t — >• t~ confirming the blowup at t c for all values of the shell model 
parameter e. Figure [TUb gives an example of the blowup solution u n (t) for e = 0.01. 

Therefore, in the case of quasi-periodic and chaotic attractors, solutions u n (t), 6 n (t) blow 
up in finite time. The asymptotic structure of these solutions near the blowup is universal, i.e., 
they scale as (I78p - (l80p with the single universal exponent y and the dynamics of renormalized 
variables v n (t), i? n (i) is described asymptotically by the attractor of the Poincare map. 

10 Conclusion 

We considered the inviscid shell model of convective turbulence suggested in [5] and studied the 
blowup problem, i.e., formation of a finite-time singularity for smooth solutions of finite norm. 
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Figure 10: (a) Dynamics of renormalizes variables v n [r) represents a chaotic wave moving with 
average speed a toward large shells n. (b) The corresponding solution for original shell variables 
u n {t) blows up in finite time t c . The blowup is hard to recognize visually. However, it satisfies 
the scaling laws ([7H|) - ([80]) with the universal exponent y, see Fig. [9b. 

The blowup condition similar to the Beale-Kato-Majda theorem is derived. The constructive 
method for analysis of blowup is developed, which exploits the Poincare map introduced for a 
renormalized system. This method identifies the blowup and its internal structure by studying 
attractors of the Poincare map. 

It is shown that solutions of the shell model blow up in finite time for all values of the param- 
eter, which controls the process of energy transfer. Depending on the value of this parameter, 
fixed-point, periodic, quasi-periodic and chaotic attractors of the Poincare map were detected, 
and the corresponding structure of the blowup was described. In each case the asymptotic form 
of blowup turns out to be universal (independent of initial conditions) and described by the 
specific scaling law. 

A similar renormalization procedure can be used for continuous hydrodynamic models. For 
the inviscid Burgers equation, an asymptotic solution for renormalized Fourier variables repre- 
sents a traveling wave [29] and describes the universal self-similar form of blowup [16]. This 
construction is analogous to the traveling wave corresponding to a fixed-point attractor, which 
determines self-similar blowup structures in this paper. Our analysis demonstrates that the 
blowup structure can be quite sophisticated and hard to recognize, e.g., the high-period or 
chaotic blowup structure. Thus, using special renormalization techniques is crucial. This knowl- 
edge may be useful for studying the open problems of blowup in incompressible Euler equations 
and fully inviscid Boussinesq equations. 

The obtained results can also be useful for explaining spectra of developed turbulence in the 
inertial range. It was shown in [28] that the blowup phenomenon is the only source for anomaly of 
turbulent spectra (deviations from scaling exponents given by the phenomenological Kolmogorov 
theory [19]) in the Sabra shell model. In this case self-similar blowup structures develop in the 
inertial range and dissipate in the viscous range, leading to correlated phenomena called the 
instantons. For the shell model considered in the present paper, analysis of instanton statistics 
in turbulent dynamics is a difficult problem due to nontrivial, e.g., chaotic blowup structure. In 
this case, identification of the instantons must rely on the knowledge of the blowup structure. 

Since scaling properties of the blowup is the important issue in the above applications, we 
summarize that there are two types of power scaling laws describing the blowup in the shell 
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model of convective turbulence. The first type corresponds to the interval —1.255 < e < 0.015 
(periodic, quasi-periodic and chaotic attractors) and described by the power laws 

Un~A*-\ 0n~h%-\ t~k-y (81) 

with a single exponent y. Velocity and temperature fields are essentially coupled in this case. 
For y = 2/5 expressions (181 h reduce to the Bolgiano-Obukhov scaling law based on dimensional 
analysis. Since y ^ 2/5 for most values of e, see Fig. |5k, the blowup scaling is anomalous. 

The second type corresponds to the intervals e < —1.255 and e > 0.015 (fixed-point attrac- 
tors) and described by the power laws 

u n ~ki~\ e n ~kly- z -\ t~k-y (82) 

with two independent exponents y and z. The velocity and temperature fields uncouple in 
this case: the scaling exponent y is determined by the velocity field only, and the exponent z 
describes small temperature perturbations. 
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